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ABSTRACT 



The analysis of scattered electromagnetic radiation is examined in 
terms of I'ie theory establishing the detailed form for the scattering 
functions. These parameters are used to define the criterion i/hereby 
the composition and number densities of an atmospheric aerosol population 
may be assessed. Laser radar, or LIDAR, is the proposed remote sensing 
device and the governing system equations are developed. The problsm of 
data inversion is surveyed with emphasis on smoothing methods, statistical 
analyses, and iterative teclmiques. A discussion of the numerical 
stability of the solution is also presented. On the basis of the Mie 
model, inversion of the data provided by the LIDilR probe is given as a 
rationale to effect an estimate of the true oarticle number distribution 



function in an atmospheric cell. 
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I. riTRCDUGTIOI-y 



The advent of sophisticated optical devices in both the scientific 
and defense communities has pointedly demonstrated the need to better 
understand the primary propagating and interacting medium characteristics 
of the earth* s atmosphere. The intent of this paper is not to deal with 
the purpose of such devices, whether they are atmospheric pollutant mon- 
itors or infrared threat detectors. Rather the concern hers is the use 
of these devices in the real world environment and the limitations im- 
posed on them as a result of atmospheric optical signal degradation. 

This degradation is the result of manifold scattering effects such as 
Brillouin, Rayleigh, Mie, fluorescent and Raman scattering and includes 
the atmospheric constituents* differential absorption of optical 
radiation. CCl 

The atmospheric scattering and absorption centers include the air 
molecules, suspended particulates (aerosols) and water droplets in fog, 
rain or hail. The presence of aerosols strongly determines the trans- 
mission characteristics of the atmosphere and may even dominate the 
propagation characteristics# See Table I and figure 1 for a comparison 
of particle sizes and the range of number densities. 

At the small end of the size spectrum are the air molecules them- 
selves with radii of 10 ^ra and a concentration of 10*^^ per cm^. The 

large end includes rain drops with radii of 10 to lO^^m and a concen- 
-5 -2 3 

tration of 10 to 10 per cm , The **small** particles, defined by a 
radius much less than the v/avelength of radiation used as the probe, 
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Number density of particles dll/dlog r (cm 



Average size distribution of nuclei in 
continental edr masses 




Figure 1 



c 



obey the Rayleigh scattering model, exhibiting a cross section that 
decreases proportionally as the particle radius to wavelength ratio to 
the fourth power. It is this fourth power law for molecular scattering 
that Rayleigh used to explain the blue color of the daytime sky. 



TABLE I 



Particle t’/nes 


Particle Radius (/<m) 


Concentration (cm 


Aitken nuclei 


-3 -2 

10 ^ to 10 


10^ to 10^ 


haze particles 


10“^ to 1 


10 to 10^ 


fog droplets 


1 to 10 


10 to 10^ 


cloud droplets 


1 to 10 


10 to 300 



initial nuclei present in the atmosphere from the earth's surface 
and/or meteoric sources. 

For ’.javelengths larger than those of visible light, this scattering is 
very small. As the particle size approaches that of the radiation wave- 
length, the complete Mie theory is required to explain the scattering 
characteristics. The dependence on v/avelength is highly oscillatory and 
very complex. As the particle grows still larger, its behavior in the 
radiation field approximates that of a large solid object casting a 
shadow field in proportion to its geometrical cross-section. The Mie 
theory includes both Rayleigh scattering for small radii and the "non- 
selective" scattering in the limit of large particle radii. The theory 
is covered in some detail in the next section of this paper. 

Once one has both the atmospheric density and data on the aerosol 
composition and size distribution as functions of altitude, location, 
and time, the transmittance can be computed on a real-time basis by some 
rather sophisticated computer models. 

Aerosols are colloidal particles dispersed and suspended or falling 
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slowly in the gaseous mixture of the atmosphere* Some examples are 
smokes, haze, clouds, small droplet (less than 1/«m) fogs, and finely 
divided soils* These aerosols have as many sources as they have diff- 
erent forms and compositions* The chemical reaction of various atmo- 
spheric gasses whose minor components are nitrous oxides, terpenes, and 
hydrocarbons can produce solid particles* The presence of combustion 
products from inefficient burning, which may be locally critical, and, 
on a much larger scale, from forest fires and volcanic action serves 
only to increase the atmospheric load of particles that may take years 
to settle* Over large bodies of water, the dispersion of solutions into 
the near atmosphere region brings about the formation of fogs and rain 
as well as clouds* As a class, aerosols may either help or hinder the 
researcher* These particles may be used as tracers enhancing a received 
signal to detect atmospheric motions, distributions, and chemical com- 
position* However, they primarily provide a large, random source of 
environmental noise* In this latter mode, the particle’s presence tends 
to confuse the received signal at the sensor and provide noise that 
originates outside the system electronics, maJcing it difficult, if not 
impossible, to supress* ^6} Historically, aerosols have been divided 
into classes to • facilitate their theoretical modeling* These models 
involve several free parameters that ere designed to "best fit" the 
observed data* The classes are: 

1* Maritime - This group is characterized primarily of sea salt 
\vdth some water of hydration and a clilorine concentration slightly higher 
than sea water* Salt particle concentration decreases rapidly above a 
height of about 500 m* 
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2. Rural - The largest number of components are silicon, iron 
and some sulfates with approximately 10^ to 15% organics* Of this total, 
about l/3 of the aerosol is hygroscopic, whereby air moisture is readily 
absorbed. This same class is also found at sea at higher altitudes vdth 
a laJTge particle density falling off rapidly as distance at sea 
increases. 

3. Urban - Combustion products and the issue from industrial 
processes maize up the bulk of this class. 

4. Troposphere - This is an atmosphere division lying above the 
surface boundary layer and consists of the iraral/maritime classes with- 
out large particles. 

5. Stratosphere - Primarily (30^ to 90f^) consisting of the large 
sulfate ions and particulates, this group also shows a periodic large 
influx during volcanic activity. 

The typical atmospheric analysis technique is to use two-ended sys- 
tems requiring a source of radiation of Icnovn properties, a propagating 
path tlirough the atmosphere of specific length and reasonably loiovn 
conditions, and a stable receiver at the pathos end. Exp'eriments such 
as these are able to measure aerosol size and composition distributions 
to check out the various composition and scattering models. They also 
predict the performance characteristics of a particular optical device. 
However, two v;ell- separated platforms are required v;ith the attending 
problems of phasing and communications. Also, only low altitude, hori- 
zontal v;ork can be done. A one-ended system would eliminate the need 
for one of the platforms and enable a full altitude/azimuthal capability 
in a real-time environment. There is a system that can easily be 
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adapted to the one-ended process and is available today. The laser 
radar, LIDAR, system is similar in operation to Radar with optical 
wavelengths instead of radio waves. In fact, the acronym is somewhat 
the same, meaning "Light Detection and Ranging". Conceptually, a laser 
pulse is emitted into the medium to be analyzed. This pulse interacts 
v/ith the components of the medium and is both scattered and absorbed in 
distinct ways depending on the material type and density. Some portion 
of this pulse is returned as an echo to the point of origin where a 
telescope-lilce receiver is coupled to a photodetector that converts the 
received optical radiation into a proportional electrical signal. The 
pulse return timing is gated to the range of the atmospheric cell under- 
going analysis. The intensity of the echo, the echo’s polarization 
change, and any measured frequency shifts all undergo involved processing. 
The outputs can consist of composition, type, and size distributions. 

Some state of the art applications of the technique include; metero- 
logical LIDAR investigations; middle atmospheric LIDAR studies; tropo- 
spheric chemistry and diffusion research; studies of atmospheric 
propagation and radiative transfer; absorption, Raman and fluorescent 
spectroscopic applications; and pressure-temperature measurements of 
atmospheric layers. A recent theoretical paper indicated that a 
series of measurements involving optical scattering and absorption can 
yield about 7 or 8 independent pieces of information concerning the 
aerosol size distribution within the analyzed cell, covering a range of 
particles with radii from 0. 1^m to 5/<m. f2(3 

Consider such a single-ended LID;\R analysis system on board the type 
of mobile platform in use today by both scientific researchers and the 
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military. The optimum performance of a specific optical scanning/de- 
tection/ tracking device could be evaluated in a real-time scenario with- 
out the need to resort to generalized models applicable only in 
standardized areas. 



u 



II. GSIISRAL SCATTERING TIEORY 



Scattering is the process whereby an incident electromagnetic wave 
loses energy and some portion of this energy undergoes re-emission. In 
general, the scattered wavefront has a different propagation direction 
than the incident wave. See figure 2. 

For particles with a low number of electrons, the scattering can be 
computed using quantum mechanics. In this sense, the incident v;ave can 
only interact with the stationary states of the scattering particle's 
electronic system within a small bandt-;idth of frequencies about the 
allowed transition levels. One can discuss Rayleigh scattering as the 
absorption of the correct energy photon and the subsequent re-emission 
of a photon of essentially the same energy. An energy level diagram 
would look like figure 3. 

Likewise, Raman scattering involves a definite frequency shift be- 
tween the incident and scattered radiation. If the scattered frequency 
is less than that of the incident, the process produces the Stokes 
lines spectroscopically. Aiid when the reverse is true, the anti-Stokes 
condition exists. 0s2 See figure 4* 

However, for molecular and larger particle scattering, the collection 
of electronic states becomes great enough to let us deal with the prob- 
lem classically. The particle is considered to be a solid dielectric 
with a complex index of refraction. 

One of the most widely studied problems in diffraction theory is 
electromagnetic and acoustic wave scattering by a sphere. The problem 
of diffraction of a plane-monochromatic wave by a homogeneous sphere 
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characterized by its radius and complex index of refraction was rigor- 
ously solved by G. Mie in 190S, His analysis involved the solving of 
Maxwell's equations for the electromagnetic field when an incident plane 
wave interacts with a material interface. '.Then the discrete boundary 
separates media of different optical characteristics, a scattered wave 
is generated. 

Today, an extensive bibliography of approaches to this problem exists. 
The author will follow the analysis presented by Kerker (1969) 023 
throughout this section in order to present a consistent derivation. 

First, the salient features of the Maxwell theory will be listed, lead- 
ing to the wave equation. The solution satisfying the appropriate 
boundary conditions will be obtained; the resulting infinite series 
solution not only completely defines the scattered wave, but gives the 
electromagnetic states of the parti cules' interior as v;ell. The single 
particle scattering results 'will be compiled and extended to cover 
collective scattering. 

Max’well's equations, in rationalised MK3 units: 



Coulomb' s la’w 


V-Q 


= 


(1) 


jimp ere 's law 


rx w 




(2) 


Faraday' s lav; 


Vx £ 




(3) 


absence of free 


magnetic poles 


V ro 


(4) 


The sjTnbol denotes 


a vector quantity. 


Folloving are 


definitions 



of vectors used in derivation: 

- electric field intensity (volts/meter) 

^ - magnetic field intensity (amps/meter) 

2 

D - electric flux density (dielectric displacement, coulomb/meter ) 
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..®t V 




2 

^ - magnetic flux density (magnetic induction, weber/neter ) 

2 

J - electric current density (amps/meter ) 

- electric charge density ( coulomb/ meter‘d ) 
t - time (seconds) 

These equations uill define the electromagnetic field config’Jration 
for all points in spa.ce and uniquely determine the field vectors giver. 

(5) 

( 6 ) 

(7) 

are: 

jjnductive capacity, farads/meter) 
inductive capacity, henrys/meter) 
conductance, 1 /ohm-meter) 

e= 3.8542 X 10“'^ farad/m 
10“’^ henry/m 

As a consequence of the invariance of the Maxwell equations, the quantity 
has the units of velocity and, again, in free space: 

C= X 10® m/sec ( 6 ) 

This last quantity is the speed of light in vacuum, ’./hen the quantities 
5,^, andc^are independent of direction, the region in space they 
characterize is said to be isotropic. Two core equations ere required 
to complete the picture. 

Lorents Force £ = X §) (9) 

I'e’,rton's second law of motion ^ (13) 



the media characteristics belou, 

D = € E 
&= 

^ E. 

The media parameters thus defined 
& - permittivity (electric 
- permeability (magnetic 
^ - conductivity (specific 
In free space (vacuo): 
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'fnere: ^ - force on charged particle (newtons) 

- velocity of charged particle (ui/s) 

31 - mass of charged narticle (kg) 

2 

a - acceleration of charged particle (m/s ) 

/>/ 

All macroscopic, linear, electromagnetic phenomenon are described by 
equations (1) tlirough (10) vrlth the associated quantity definitions. 

The clas33.',cal d^^raxiics of interacting charged particles and electro- 
magnetic fields are thus specified. 

A certain degree of convenience is obtained using the definitions of 
the tvro following auxiliary vectors: 





(n) 




02) 



Tnere: ^ - electric polarization (c/m ) 

^ - magnetic polarization (V^) 

One can see that, in free space, these quantities vanish identically and 
they may then be considered as a measure of the effect of natter on the 
locai fields. Some dimensionless quantities may also be defined: 

- §/e, (13) 

'."aere; \(^ - specific 6 (reia.tive dielectric constant) 

- s?ecificy4<(relative magnetic permeability) 

The media may further be described with: 

05 ) 

0 « 

I'/here: electric suscept ability 

magnetic susceptability 

The previous set of equations can be combined to obtain: 
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Xe = Ke -I (’7) 

Xm= K'rr,-! 08) 

Maxuell ' s equations apply to regions v/here the aedia may be repre- 
sented by the continuous functions and O', ’.fnen a medium discon- 

tinuity is encountered, such as the conditions across a regional 
interface, the fields must satisfy a specific set of relations at the 
boundary. See figure 5. Let; 

2 

^5 - surface current density on interface S (A/^^ ) 

- surface charge density on interface S {c/a) 

The Maxuell equations (l) - (4) are cast in point form for each 
differential point in space. One can apply various integral theorems 
to change the same equations to integral form and use them to define the 
conditions that the fields must satisfy by crossing the interface. [^1'^ 
The boundary conditions obtained are inherent in the Maxvell equations. 

Boundary conditions: 

(a) the tangential component of ^ is continuous 

(£2 ~ ^ £ =0 (19) 



(b) The tangential component of K undergoes a discontinuity equal to 

a= Is (20) 

(c) The normal component of D, undergoes a discontinuity equal to 

2 = ^5 (2’) 

(d) the normal component of ^ is continuous 

(g2-g.)-R, = O (22) 

For all cases of dielectric scattering, it is assumed that there is no 
free surface charge or current density. This assumption begins to brealc 
down for large particles in the presence of storm activity. However, ;;e 
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interface 3 



Figure 5 
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’.■all, for simplicity, assume that ^ = 0 and(^^ = 0. Therefore, the 
tangential components of E_ and H and the normal components of D and ^ 
are continuous across the boundary. 

Charge conservation is included in the Max’.jell theory and is written 



as: 

continuity equation (23) 

VJhere: - volume current density (Vm ) 

- volume charge density (c/m ) 

Throughout the follo’wing derivation, a feu vector identities will be 
required. 



ITote: continuously differentiable vector field 

0 - continuously differentiable scalar field 

V • (V X A ) = o 

V = vV 

V = o 

V-{0^) = A • (V0) + ^2^(v-A) 

V X (vxA.') = vCv-d) - v^/i 

If the fields described by (l) - (4) are invariant under a gauge trans- 
formation, then they are v;ell defined by the preceding formalism. 

VJe define: 




vector potential field 
0 - scalar potential field 
in the following relations. 

B : V XA (25) 

Compare this to equations (4) and (24). Using equation (25) in (3), we 
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get the following: 

How, refering to (24), 



=> E +• 






■i= -VJZ^" 



, thus 



E - 






■4-V^ (26) 

A and 0 are not unique. Let g be any scalar function, then the new 
potential can be defined as: 

A ' = ^ - VS (27) 

0'r $2^ + (28) 

Now the fields E and 3, are recomputed, 

B'=rx (d-^3) = B- Vx(r9')'B 
£'=‘itC4-va)-v(^+-|r3)= '^A-'!'sf+ir(’’3)-r 






The physical observables are the fields ^ and ^ and it can be seen 
that the gauge transformation defined by (27) and (23) leaves these fields 
unchanged. Thus, the Maxwell equations are gauge invariant, 

'.7e now turn from formalism to a more practical aspect of electro- 
magnetic theory and introduce the Wave Equation: 
r X ( 6 ): VX (r X E) i ^ (V X B) 

Note that the time and space operations of differentiation are commutative 
for these continuously differentiable functions, 

(0, ( 2 w): r (v-£ ) - r^£ = '-fr 

The assumption is always made that all space and time derivations of the 
parameters 6,^, anddare either slowly varying \/ith respect to the rate 
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of change of the fields themselves oc identically equal to zero. 

=■> v(vi]-v^i= 

CO, ( 2 ), ( 5 ), ( 7 ) ; V. f&£) = e V. E = 

E = /y£. 



=■> V 



r>rxR= cr B + ^ 



=^> V C/^/£ ) - V ^ g = y^ig ^ 

The same argument also applies to the quantity /^/8 concerning its rate 

of change with the rate of change of the fields. 

^ i, yU.21 E 

-kt ^ £Ti- ^ 



— > - kt 

2£ _ - £'/'21 E ^ o 



Therefore: - O (29) 

This eo^uation, (29), is the uell-kno’.jn nonhomo gen eous equation for damped 



wave notion - the '.lave Equation. SLmilarly, V X (2); 

vx (vxH) = 

V • (y^ ^ \i - O i=> V • H = O 

(3): -V"!i = 

C^) : - V ^ H 



rr 2- , / - _ W H 



at' 






— > V^ii- = O 

The relations for non-dissipative media are obtained by setting o' - 0, 
Both equations (29) and (30) actually represent 3 equations each. Each 
component of the vector field must satisfy these equations for the vector 
field itself to do so. Let the symbol u correspond to any of the six 
components of the two fields E, and The scalar wave equation can be 



expressed by: 

V - cf /^ ^ - ^/^djL = Q 






(31) 
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The fields described by (31 ) represent a transfer of energy in both space 
and tine. The electromagnetic wave has a time periodicity that sug- 
gests a fora for the energy flux crossing a unit area per unit time: 

S = E H (32) 

Equation (32) is the Poynting vector equation, whose tine averaged 
magnitude over the period of the wave yields the wave intensity, 

I ^ '/r dt (33) 

whereTis the period of the wave. The field symbols E and H in (32) are 
the real parts of the complex fields v/hose components satisfy (31). 

The form of the scalar vave equation suggests a harmonic time de- 
pendence of the form: 

U (iC.'t) = (34) 

\jhere o) is the angular frequency, i =V(-1 ), and r is the usual three- 
dimensional position vector. Note that: 






LujU 



_ Z . I 

and = u) U. 



d-t 

then, (31): V% - cco U ^ co ^ U O 

=> - L/^ cr uj) u = O 

Define the propagation constant k by: 



(35) 

then: V U = O 

Equation (36) is the Helmholtz Wave Equation, teother form of (35) is: 



K = e<j — t X3| 






'4 

V2 



(37) 

(38) 

(39) 



It is the/3, term that determines the damping undergone by the wave per 
unit path length; the energy is dissipated as Joule heat in the medium 
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v/henCS' >0, 



Recalling that the quantity 



V - 



( 40 ) 



is a velocity which, in free space is c, the relation for an electro- 
magnetic wave in a nediun described by light speed V is 

^ (41) 

where: X - wavelength 
■0 - frequency 
V - velocity of the wave 

Also, 2Tr\) = u> (42) 

and this is the angular frequency, u> , mentioned earlier. In free space, 

(43) 

The index of refraction, n, is defined as 

(-U) 



=c 



n = c/v 
,): n 

=> n X = Aq 



(41), (43), (44): n= 



(45) 



It is clear that A.© is the wavelength when the v;ave is propagated in 
free space. In the sane medium. O' = 0 implies that 

/3, =1 0 => '-''I = to - u>/c 



=> Ko= ZV-A 



(46) 



This is the free space propagation constant. 

The complex refractive index, m, can be defined as 



m 



= k/Ko= n (I - iK) 



(47) 



where: hC - index of absorption (extinction) 

As an example, suppose that there is a plane v;ave propagating along the 
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z - ajcis whose x - component 

LUJ^ 



( 47 ) : ^ 

—V £;< =■ A e 
■=>£x A e 



- ^-cr\f<!e>0 - c. X) Z 

iwi .ink^z. -nk^Kz. 

& 6 ° 

i((vi- nKoZ) ^-nK, Xz 



The attenuation of the v;ave is seen e;’:plicitly as e 



-rf' X: 



■in e:coansion 



of this enample shows that the resultant plane wave solution leads 
directly to the Stokes paraneters. Assune the propagation direction is 
s, as stated earlier, and that the nonoclironatic wave is traveling 
unbound in a source free isotrooic nediun. It is assumed that the wave 



by completely specified in 


the ::y-plane 


'Wave. Thus: 






Aik . AJk 

5 X ^ y 



(34), (36) 



= O 

li^t 

e =0 






Substituting a trial solution, f - ae^", 

=> ak^e‘’^+ K^ae = o => U’‘= b = i (k 

=•> f(z) = 

From the definition of k, -Hcz is selected as the vaJ.ue. 

^ 4.'^ - '■ 

=y U 12 :, t; = ae 

This equation describes the fora for each of the six vector components 



of E 


and H. 


Recalling (l 






d£<i 




0 ) 


dx ^ 




3Z. 


But, 


ax = 


0 ) 4 — 


0 '> 




£ 


i {(uji’ - 



= o 



and U ^ £ 

= ay e 
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mm* 



( 3 ), ( 6 ) 



V xe: = -> _ liiy 

at ^ cPt az'^’ 



y - 

hi ^ - <5z 



-^hH 

"> 6t 



=> 'Tt' V X *> -TT- Z^O => (^ 2-0 

It is quite sufficient to describe the behavior cf ^ as the value of ^ 

can be readily obtained using equations (l) through (4). Given the form 
of u(s, t) : 

SH 



h£ 

=> 



=> i U> rfy 
:.OT,e: ^xy'hl^J - <-tj/ K 

The term/^^/k can be reduced using (35) to the folloid.ng: 





3 


My 


z i a) 


II 

1 


3 


-~y 

a z ^ 


= -iA 


•t i Pc £if 


■zzy 


11 

X 

Tt 


^ 1 


~ ik£^ 








^ a; ^ 



K " y^uj \ ^ to ^ 



Defining the intrinsic impedance of the medium for plane naves as 



■ (Ar) 






(43) 






-E 



-V 






Hy = 

finally, 3^ for the plane wave is: 

t. (a)t - Kz ) 



E^= Ae 

£,= 0 



i (ct)± - Kz. +c^- ■) 



(49) 



where and ofiare arbitrary phase angles to account for the general his- 



tory of the wave. The !I i;ave is also determined by equation (49): 

t _ i(u)t-/<z foTj) 

^ X “ ~To ^ ® 

l Cidi -^2, -/-J'l ') 






(50) 



= O 

;/e see that: 






£y £ 
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( 51 ) 



=> E • H = O 

Thus, T and "" are *oernendicular to each other for the olane wave 

solution. 

The concept of polarized light can novj be introduced with the help 
of equation (49),^",, and^. The general plane wave solution describes 
an elliptical rotation of the 1 vector in the olane pemendictilax to the 
direction of propagation. Defining the phase difference as«/ : 



^ ^ Sy - S, 



(52) 



The representation of the state of a polarized wave can then be uniquely 
deternined using the Ttokes parameters, which are: 

bJ (53) 

s, = = So cos 2 v* cos 2 X (54) 

2 f\b C 05 if = Sq sm z 'f' Cos z X (55) 

Ss= 2A6 S1 NcT= So sin Z X (56) 

The lengths and 3, and the angles and X are illustrated in figure 6. 

The quantity (cot - kz) can be eliminated from (49) using equation 
(52), where E and refer to the real parts of the complex v/ave: 

This is the equation of an ellipse. The ellipse in figure 6 is drawn 
"right-handed", that is, the head of the electric field vector folio v;s 
the arrows when viewing the wavefront face on. In order to exhibit right- 
handed polarization, sinS>0 and 0^X4fl"/A* Linear polarization exists 
when: 

S = j= 0,11, i 2, t 3, - • (58) 

and circ^.ilar polarization corresponds to 

fej-r, j=tl,±3,ts, ±7, ... ( 59 ) 
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Pictorial definition of the 
quantities ^ 




2A 



tan X = i V«, C%^ X $ 7v), (0 4 1T) 



Pigure 6 
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The Sookes paraneters satisfy: 

S© = ^ I S 2 + ^3 (^0) 

and a complete state of polarisation is obtained from only three in- 
dependent parameters. 

Natural light is defined as the state where the intensity of light 
is the same in any direction perpendicular to the propagation vector, lu 
The phases of the various components in natural light are randomly 
varying so that no correlation exists between them. 



=r> 



n _ 



= Sj= =0 



' , - ~ - w ( 6 l ) 

The combination of eUiptically polarized beams can be accomplished using; 



Vi ) S, - ^ S,i. J S2- > “^3 ‘ f- ^31 



( 62 ) 



Then for non- identical polarization states for all beams: 

sj > s? S| ^ S 3 ( 63 ) 

Partially polarised light can be described as a superposition of natural 
light and a set of incoherent beams: 



s - s 

J I 



S, = sf , Sj -- 5 



P 
Z ) 



and 






(64) 



Finally, the fraction of polarization, F, can be defined as the intensity 
ratio of the polarized part of the wave to the total wave superposition: 

^ ' ^0/(60 ) (o^F^l) (65) 

It can be shown that the parameters of a scattered wave are linearly 
dependent on the parameters of the incident wave. This relation 



is written as follows: 

- 



I, = (l/RhSIi 



( 66 ) 



where: R - observation distance to scatterer 

i - incident wave 
s - scattered wave 
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The nodified Stokes paraneters are defined first by writing equation 



(49) as; 



’.Oit 






LUit 



tnen. 



e, =e,2 

I. = 1E,1 

X. -Uz 
U= 2 Re 
V z 2 [E| ^2*} 

The 4X4 Stokes matrix ^ looks like this: 



(67) 

(63) 

(69) 

(70) 

(71) 



a - 



2 Ir»nCf|(f2|) 



V_ 



A 

-I-CfAt-fJzO 

2 Irr.Cf.Jjt) X,vft„fzt*fafzO Re tet-hfzt)/ 



f„r 


Ifj' 








2 Re(f,i4?) 



(72) 



where the f^^. are defined in terms of the scattering functions given in 
&0- For spherical particles, the quantities fj2 = ^21 “ ^ ^ 

contains four constants. 

Having laid the foundation for the Kie treatment of scattering 
phenomenon, let us now consider a plane electromagnetic wave incident on 
an isotropic, homogeneous sphere of arbitrary size. The object interferes 
;d.th the free propagation of the wave and secondary waves are generated 
that are predictable by Maxwell's equations and the boundary conditions; 
and these secondary v;aves form a scattered radiation field. The incident 
wave should be linearly polarized for two reasons. First, if the light 
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is unpolarised, then it can be resolved into two linearly polarised 
conponents each of which act independently of the other. One only needs 
to consider the interaction of the sphere with a plane (linear) polarised 
wave and then apply the principle of superposition to obtain the total 
effect. The LIDAR laser pulse is plane polarised, by definition, and 
this is the wave-type that is considered for tropospheric probing. A 
geometric diagram of the problem is shown in figure 7, suggesting a 
spherical coordinate system as the most convenient for analysis. 

The fields of electromagnetic energy in the region of space filled 
by the particle and the incident wave are resolved into three conponents: 

a. the incident wave - E., H. 

b. the narticle interior wave - E , H 

-r^ >rr 

c. the scattered wave - E , H 

For these quantities to properly define the fields, they must each 

satisfy Manwell’s equations (l) - (4), the IJave equations (29), (30), 

and the boundary conditions (19) - (22). The internal field E , H 

must natch the external field E. + E , H. + H in accordance with 

^ *^1 -^s 

these relations. 

KerkerQ^et al. , suggests that the derivation deal vdth the scalar 
wave equation, (31). Two new functions may be introduced, the electric 
Hertz vector,)^, and the magnetic Hertz vector, Tl^, and are defined as 



follows : 




(73) 

(74) 

(75) 

(76) 
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Scattering of 



a plane wave by a sphere 



r7 

Lx 







incident plane v/ave 



Fig’ire 7 
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Suitable transfornations may be applied to define other Hertz vectors 



as long as the fields are left invariant. The vectors thus defined must 
satisfy the vector wave equation: 



dipoles. Therefore, two vectors can be used to describe all the fields 
present. 

The transverse magnetic ’.;ave (TM), also called the electric wave, 
has a zero value for the radial component of its magnetic field intensity, 

- 0, Likewise, the transverse electric wave (TE), the magnetic wave, 
has a zero radial component of electric field intensity, = 0, The 
TM wave is conceived as being the result of oscillating electric di- 
poles in a spatial region, and the TH v/ave as the result of oscillating 
magnetic dipoles. The total solution is a construct of these tv;o 
wavefronts. 

The Kertz-Debye potentials are: 



These potentials are the solutions of the scalar wave equation. Solving 
(79) and (SO) for the potentials and applying the appropriate boundary 
conditions, one can then find the field vectors describing the TM and 
TS waves. The total field may be derived by the direct addition of the 
component waves. 

The follo’wing equations for the total field vector components can be 
obtedned by applying these concepts to the above relations. In spherical 




(77) 

(78) 



where the^ and result from distributions of electric and magnetic 




(79) 

(80) 
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coordinates, we have: 

Er -- E|,^£2r = ^ Cr%) + K^CrTi:) + 0 

ar3© ^'’"‘^■'■l^z'rsin© C’'"z') 

Eji = ^2^ - r Sin © JrT^ £''^1'^ “ ^2'f © C'* '*^1^ 

Hr - H,r +H2r " ^ ('’‘n'i)+ K^CrTt^) 

riin© ■^p' (r + r J^r3e 



H;:^rHlgf + H2/= K, -j: -|^(rTr>Vi:4;7r© 

vdth the propagation constant defined in the follo\ri.ng fore by (35): 
1:^ = "K, K^{c7), '.;hereK|=- LoJC' + o' (SS), K 2 = ^ dJ (c9) 

has been dropped as all nedia are assuned to be non-nagnetic. 
Recalling the scalar wave ecuation for a function vrith harnonic 



(31) 

(32) 

(83) 

(84) 
(25) 
(36) 



u . j , iwt 

iine depenoence, e , 

iujt 



we nave: 



(36) + K^a = o 



(34) U = £(!)£“ 



The spatiol part f then sewtisxies: 

2 



K -PCc) - O ( 90 ) 

The Hertz-Debye potentials are constracted such that they are solutions 
to tliis ;:ave ecnation. Let IT (r, 0 ) " f(r) represent oither 

potential, 7T, ov\: 

V^TT -h k^iT o (91) 

This equation is the one that must be solved; in sphericsl ccordina.tes : 

\ f 1 ^ c)1T\ J </_ z . . 

r C^T*’ ^ *0 r*sm© de ( 6 Q/^ r^3ir\"^© K 'TT = O 

It ca:i be solved using the technique of separation of variables, whore 
fTis assuned to be a product function: 

Tr(c) = R(<) ®Ce) (93) 



-> C'"^) = ® ^ r R 



(qr\ 
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=> ® '°°' 
(92! beoonss;^ ^ I5 (siM S |a) 4.-^^ 



■Tv'll -Ki 



2 

■ti‘:i:/ bv r /iT : 



'ir.GG 'll" j then '.;q 

-L ^ i_ Ain e . 

R 3r^ ® SIN e ^ i+ 

r 



r 



S I . ^ 2 

~h~^ - “K r 

K+KV^ P A/swea^J Hi~I 

r -'g)siN© a©V d©7'*isnj^e ^ 0^_j 



rc.n) 



rRCr) = O 



R ari 

Tl'_Is equation is equal to a. constant which nay be cheched by tal:ing the 
derivatives of both sides with respect to r, © , {ZJ . Tliis constajit has 
the forn of n(n + 1 ) : 

r R + = nfn + O 

Thus the radial equation is: 

RCr) -f 

d 

Using (9'f), equation (97) becones: 

2i aVsjv©a(S>'i j d'^I. ,h 

®SIN©^( 7ey + ^:s/N*© = -r\(n + i) 

rsiwe i ® 9® \ , „ / .,^ «5,w _J 

LTSr^ I©( n(n^i)s.M 0J= ^ 

By the sane differentiation aTgu-ment, each side of this equation is 

2 

equal to a constant, m : 



ua^ion is: 



(98) 



Thus, the Q -angle eouation is: _ 



(99) 



Using (98) and (99), (97) talces the forn of: 

1 T¥i- - ^ 

Thus, the 0 -angular equation is: 



(100) 



In the above equations, n nay talae on any integer value and n nay be 
any one of the follovring set: -n, ..., 0, ,.., n. Equations (99) and 
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(100) conbir.e to forn the Spherical Hamonics vrhich.are defined: 

C {®,^) = fln(E)n (®) 



n < 

/here A is a nornalization constant to ensure orthonornality of the 



n 



function set of Y. 

The solution to (lOO) is a straightforuand coabination of siiies and 

2 . 



cosines as a is al’,'a.ys positive: 

C 05 + bm 5in (m 



(102) 



liquation. (99) is the fauili-ar Associative Legendre equation having 
the re.^-uLar solut: 

(103) 



©ce) = 



(rn) 

where the 7 )^ (cos © ) are the associative Legendre ooL,;Tio:niaL3. 
Using Rodrigues^ fornula : 

•n (C6S0) - (zn.) Cl-Xj ^ ^ rrtTa (lC4) 

where :: = cos @ and the negative values of a are included, -n^ra ^ n. 
When n ~ 0, then (99) is the Legendre equation whose solutions are the 
Legendre nol'^aionials: 

'd 



P„(X)= (^) ^>'-0'' 

Ycuation (IO 4 ) may be ret.’ritten as: 



(105) 

( 106 ) 



The complete spherical harmonic solution includes the normalisation 
constant and reuTites the azimuthal solution as a comple:: e::pozienticl. 
These functions are completely orthonormal as well. Thus: 

v:(e,^) 007) 

The (-1)*^ term is a phase factor defined as the Condon-Ghortly phase, and 
may, in some texts, be dropped. 

(hote: in the remainder of this paper, I './ill use ^ instead of zeta. 
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uMch is used in the literature.) 

The Radial equation, (9£), has solutions defined in terms of the 
Bessel functions: 



(lOS) 

Xn(X'') = "yf Kr A/n^'/2.CKr) ( 109 ) 

=> r = Cn CxO Xr (no) 



These are the Pdcatti-Bessel functions, t-;here J i and 
integral order Bessel and heumarui functions. By using 
definitions of the Bessel functions, one can ’..rite the 
function in the follou-ing derivative form: 

v/here s = x + iy 



il . X are the half 
n+v 

the infinite series 
Ricatti- Bessel 

(m) 

( 112 ) 

(113) 



The function: 



X„CKr) = ^/^ 



(114) 



vanishes in the linit as kr approaches infinity, and v;ill be useful in 

(2) 

later work. The function H\i(kr) is the half integral order Hankel 

n+2 

function of the second kind. 

How it can be seen that the solution to the scalar wave equation is 
characterised by the numbers n, m: 

nrr'= rR„(r)©<”>(e) 

and (102), (103), and (llO): 

=> rtr^"''’(r,e»= (115) 



Tile general solution nay now be obtained as a linear superposition 
hese particular solutions. 

rTTCO =11 rlr<”\r,e,«<) 0l6) 



n:o mi-n 
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Equation (116) is the general solution of the scalar wave equation in 
spherical coordinates* 

The fields in and about the particle can be described in terns of 
their individual pairs of potential functions; 
a* the incident wave - rr,‘, n 

b. the interior wave - Tir, 

5 5 

c. the scattered wave - rr^,rr; 

The sphere, whose origin coincides with that of the spherical coordinate 
system and is of radius a, is isotropic and homogeneous, optically 
described by a complex propagation constant, : 

(47) ~ ^0 (117) 

The isotropic, homogeneous medium is a dielectric described by a real 
constant, 

(113) 



The relative refractive index is defined; 

K,/k^ ^ rvi, (119) 

The incident plane wave propagates along the positive z-axis and is con- 
sidered to be of unit amplitude with polarization ajcis parallel to the 



x-axis* 






( 120 ) 



In order to facilitate matching the boundary conditions for the Debye 
potentials on the surface r = a, (120) is written in the form of (Il6), 
realizing that the function Xr\(^2^^ becomes infinite at the origin 
where the field (l20) does not: 






( 121 ) 

(122) 
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The unperturbed incident wave is now well defined. Because the incident 
wave propagates along a coordinate axis (a), parameter m = 1, As a 
result, the Debye potentials all have the associated legendre poly- 
nomials of degree 1 . The internal particle fields must be of the same 
form, xd.thout ^^(k,r), and must also be finite at the origin. Also, the 
arbitrary coefficients must exist in order to properly map the conditions 
at the boundary; 



= d; ' s)c-i>s (123) 

rTT/ r Tf “i i"" dnt,CK,'-) P“(«5 s) s.« (124) 

The scattered wave must vanish at infinity without entailing an infinite 
energy and the function ^r,(kr) contributes just this property while 
maintaining the proper form for 71^ (r). Arbitrary coefficients are a- 
gain necessary: 

~ ~ ^ ^ ^r> (coi e) Cos 0 (125) 

='ih I-,''’'' 1^(126) 

The boundary conditions at the interface ensure the continuity of the 
tangential comnonents of E and H. In soherical coordinates, this trans- 
lates into : 



Eq (r = a) = (rza.) 
E^ (rzo.) — E0 C 0-) 
( r-a) = ( r -a) 




(127) 

(123) 

(129) 

(130) 



where (l) refers to the sphere interior and (2) the exterior. The mix- 
ture of terms in TTj and at first glance malies it difficult to apply 
the boundary conditions. However, by forming the appropriate linear 
combinations of (31) - (36), the boundary conditions can be recast into 
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an equivalent fora for the Debye potentials. Define: 



^ (fsiusEg) *4? 


(131) 


3 = |^(sw©He)t-^H5( 


(132) 




(133) 




(134) 


Then the boundary conditions (127) - (130) become: 

A'^'^(r=a') = A^^\r-ia) 


(135) 


B"’(r = o-> = B'*’cr = 0 


(136) 


c‘'’Cr=a)= c‘*’(r=a') 


(137) 


&"’("r--a)= £>‘^\r^a)) 


(138) 


'lote that: — TT 


(139) 


rr'^^■rr^ 


(140) 



The quantities A and C renove r and 3 and D remove rW, . The decoupled 
equations then result in boundary conditions on the Debye potentials. 
Substituting (31 ) - (36) into (l3l) - (134) yields: 

rA = ® C4' ) 

rc= f^^[)C,r7rJ (us) 

rD . f S/^J ® r/?Q r/?(3 (144) 

For equation (135) - (138) to be valid for all 0 , on the coordinate 
surface, r - a, requires the various bracketed arguments in (141) - 
(144) to be equivalent between regions (I) and (2). Thus: 





- [lcrrr.)y;i 


(145) 




- [i?crn;)ty 


(146) 


[K r<r,tL 




(147) 




= [VC.rrrd-^ 


(143) 
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using (139) 


and (140): 






--^O-cnhtb] 


(149) 


Irfirn 




(150) 




= x;‘V(E‘tir,") 


(151) 




= <'V(iThv;) 


(152) 



Equations (149) - (152) are the required boundary conditions on the 
Debye potentials. Substitute (121), (123), and (125) into (l5l): 



K 



{') 



(0 °° 



Since all terns in the series expansion are linearly independent of each 
other for each value of n and for all © , then the corresponding terns 



in the series nust be eoual: 

0) . . sx-W 



^ [r„ M- <>•" 3 '' 

y; C»<, a.} = |r„ (K^a) - a n ^ n (K 



K* yC\ 

The quaiitity ^ can be reduced by using equations (87), (89), and, 
X”i 

recalling that the frequency of an electromagnetic wave is unchanged in 



traversing the interface: 

- k!’ k7 kT Mii , 

k’ k,‘" X T K';’ ' x" ' ■ 

Applying similar argunrents to the remaining boundary condition yields a 

Ki m.Ko / 

term of the form: = — = m,/ra, =: rn 

m^,Ko ' 

The four boundary conditions then become: 

Cn Y'n 053) 

rrx [r„'(K^a)-t,^^r{(Kia.)J (154) 

Cr,VnCK,0= <’55) 

dr.'i'„(K,a) = C56) 

Bote; Y‘(nr)^ -^Y(Kr) 

This set of equations uniquely determines the coefficients a^, b^^, 
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and d^. However, only those corresponding to the scattered wave are 



inportant here. Define the following paraneters: 

=: k^a -- = 2TT (157) 

= k,o. = rnH^a, ^ 2irm,g. ^ (,53) 

A o 

where is the vacuun wavelength. 

The quantity «x is called the diaensionless size paraneter. Solve 
for c^ from (155) and use the result in (153). Also solve for d^ from 
(156) and use the value in (154). Ue obtain: 

% = (•< > Yn'C^:> -/« /"/I 

b _ 

c-^-) V'A -'fn (/3}3X C<^) 

All the field vectors are now uniquely described given the Icnovna 



(159) 

( 160 ) 



parameters ando<. This data enables the coefficients a , b , c , and 
^ n’ n’ n’ 

d to be computed and their potentials determined, Finally, the vectors 
n 

are obtained from (Si) - (36). The formal solution is complete. 
Practical light scattering observations are performed in the far field, 
or wave zone, Tliis condition is such that k^r n where n is the 
Ricatti-Bessel function’s order parameter. In tliis approximation, the 
previous relations become simpler. The as^miptotic e:r:an 3 ion of the 
Hailcel function becomes [ij : 

= (,60 
-- (162) 
The fields becone transverse due to the radial conponent decay as 
( X/r) compared to the (A/r) dependence for the non-radial components. 
Finally, using (Si) - (S6) vdth (125), (126) and (l6l), (l62) for the 
scattered wavs ortly ( the tends 'to zero for large 'c^r with re- 

spect to Jn(k^r) as does Yit'' ) yields: 
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- -L -S. 



5//V v5( C©) 






v;here the annlitude functi-ons «re: 

an+/ 



S, C®)= 7;^ [a„ 7 T„('a$ 9) tLnX, (■«!©)] 

and tho an^'ular functions are: 

ir„ c c°^ ©) = 

r„(cos e) = ^ p„^'’(nse) 

The energy flow in the scattered wave can be calculated by using 
Fo^mting^ 3 theoren: 

S ^ yz(^e>^Jl -E^He) 069) 

where - represents the complex conjugate of the quantity. The intensity 



(163) 

(164) 

(165) 

(166) 

(167) 
(16S) 



function is defined as follows: 



: (e) = IS(&)j 



(170) 



The intensity of the scattered radiation for unit incident intensity 
polarized in the ^ and© direction are: 



© ” i^rr 



cos'^ i 2 



(171) 

5 ^^ z / • 

l0 =lfp^ SIA^ pc, (172) 

2 

Note that the intensity decays as l/r as it should for this spherically 
spreading wave. The scattering plane is defined as that plane which 
contains the incident direction and the direction of the scattered wave 
(©,<20. Then i^ is the wave perpendicular to the plane and i^ is the 
wave parallel to it. The phase relation betv/een EQand S^being arbi- 
trary is the condition that defines an elliptically polarized wave. In- 
tuitively, one can think of each polarised component coming from that 
direction as being inherent in the incident wave. Therefore, the 1 9 Is 
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caused by aa incident uave of intensity cos ^ nclarized parsdlel to 
the scattering plane. 

The phase difference that exists betv;een these two parts of the 



scattering wave is: 

tan Re(Sz) +ImCS,)I,virSz) 



(173) 



IJnoolarized incident light gives rise to a scattered bean of intensity: 



= 8irV- (*-'^^0 

with degree of polarisation given by: 

i. I - L 2 



P = 



C I + I- 2 



(174) 



(175) 



Fixing the scattering plane in space to be the ys-plane, all ob- 
servations by definition tal:e place in the scattering plane. The in- 
cident electric vector has its direction in the xy-plane at an angle X 
to the y-azis. The x-a::is is designated vertical and the yz-plane 



horizontal. This is shorn in figure 8, Thus, the vertical ('/) and 
horizontal (H) components of the intensity in the scattered wave for 
unit incidence become: 





^ Xv(X) 


_x 

" WTT'' 


pt ‘i ■5"’ 


^ X 






(176) 


I® 


= rH(5f) 


- JL 










(177) 


The Stoke 


s parameters 


; are: 












So 


- ■ 4 


i, 5 /A/ 


X -h c 


; CCS 


^x) 




(173) 


S, 


- r; 


, 3/A/ 


" X - ^2 


CCS 


^x) 




(179) 


Sz 


“ V ‘z. 


3///' X 


CCS 


X cos 


cf 


(130) 


Si 


- / 


^ Z 


3//V X 


C03 


X S//V 


/ 


(1S1) 



When discussing cross sections, one is speaJcing of the amount of 
area that the particle effectively presents to the incident beam by its 
ability to spread out the incident energy. The complex index of 
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Scattering plane geometry 




Figure 8 



refraction includes absorption o.s well as scattering. Define the ex- 
tinction cross section as the sua of the cross sections for scattering 
and absorption : 

(1S2) 



By integrating the real part of the tiae averaged Poynting vector for 
the total external field: 

I = ^ (1S3) 

the total outward flow of energy is obtained. For unit incident energy, 
the various energy losses represented by the integral resulting from 



above are the desired cross sections: 



^SCQ, = -^ ^^1 ] 



(154) 

(155) 



The efficiency factor is defined as the cross section per unit 
geometrical area (ffa ) of the particle: 

Q= ( 186 ) 

Terms of the form^ X zB'-i. become since «< =ETrOL/x, Thus: 

ZfT ^ ^ ' 

Q.SCcl = ^ + (187) 

Qext = ,^:fr + [ Re (138) 

Ilote that '^sco. and are independent of the state of polarisation of 

the incident wave. 

Kerker Qij et al. , gives an excellent review of the difficulties 
involved in the numerical computation of the previously listed scattering 
functions. The functions needed to be evaluated are a , b , Tf , . 

Some discussion of the equations is given in the Appendix. 

Using approximations for small arguments, the- fourth power law 
attributed to Rayleigh scattering can be directly obtained. The first 
three scattering coefficients are: 
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a.2 



-- ^/3 C 


\. z / 


11 

- 


/4>iXlL\ ^ b- 
C 2/h’‘ +3 ) 


“ 45 i" 





(1S9) 

( 190 ) 

(191) 



plus terms in higher powers For the Rayleigh criterion«x<<l , 

only contributes. From the Appendix: 



iT, (C05 9) :: I 
T, (co5 Q) - cos 0 
(165), (i66): S,(Q) = a,TT, (cos e) 

Si(©)= i 1, r, Ccose) 

^ l<^l 1^ cos‘ & 



(192) 

(193) 



For unpolarized incident light, for example, sunlight; 

- ' ^11^/ 

(174): -1-u =aW^ H ^ 

_ q X"*- U I 
(189); Xa = Im'^t 

)\ 

1^(1 + cos^ e) (194) 

This equation accurately predicts the blue sky phenomenon as the shorter 
wavelengths scatter much more than the other visible wavelengths in 
daylight. Equation (175) yields the observed degree of polarization: 

P-. I k.r 

1 Id, 1^ + )(t,r cos 0 1 

p- I ~ cos^ 0 (195) 

I +C05 ■*• © 

A plot of this function on a polar graph indicates a maximum of 1 , 
complete polarization, at a scattering angle of 90°. This has long been 
observed in the day sky. The error associated with the troncation after 
n = 1 is 5 p if a/X is kept 0,05. This limit is the required Rayleigh 



(157): Xu 





criterion. 
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The cross sections and efficiencies represented by equations (1S4)- 
(l8S) give the total integrated differential cross section with the 
integral formed over the entire solid angle sphere about the particle. 
The differential cross section is: 



d <y c a) 

d/1 






( 196 ) 



where: - number of photons scattered into a solid angle per unit 

time per particle 

- number of photons crossing a unit area per unit time 
dn - Si-fi 0 <de ^ the differential solid angle 

The intensity is obtained by multiplying the number of photons at fre- 
quency by Planck's constant, h, as hv> is the energy per photon. 

h = X Toulfi/s (197) 

As written, equation (196) describes a scattered energy per incident 
intensity. The scattered energy becomes a scattered intensity when 
divided by the distance to the observation point, squared, for then the 
solid angle becomes the subtended area, 

d A ^ r ^ d IT- (198) 



utilizing the above, 

dl ^ 



dA 



v/here dl represents the scattered intensity over dA per 
intensity. The equations describing such an intensity, 
may be cast in the form: 

H fr 

where 0 (d) = C05 ^ or Sin ^ 

i ©) = or c ^ 



(199) 
unit incident 
(171) and (172), 



( 200 ) 
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I is the dl per unit intensity 



Thus: 






(201) 



This equation, when integrated over all angles & , ^ yields the cross 
section given in equations (184) and (l85). Hereafter, the synbol xnLll 
aean dO/d-aas defined by (201 ), This is in keeping with the notation in 
the literature. \ / \ 



( 202 ) 



This area represents the scattered energy per particle per unit intensity. 
The theory aay now be extended to multiple particle systems vdth the 
addition of the following definitions. Let n be the total number of 
particles in a region of differentiable volume, cross-sectional area and 
extent given by dV, dA, and dx respectively. 

Then : 





dV = dxdA 


(203) 


and 


n = NdV = HdxdA 


(204) 


where H - particle number density 






Assume the scattering to be incoherent among the particles. 


Then 


the total O' for all the particles is 


just 








(205) 



where ^ - scattered energy per unit intensity for one particle 
The amount of intensity scattered into dA is the scattered energy 
divided by dA. 



cf _ scattered intensity (206) 

“3^ " unit incident intensity 



.'nother useful parameter is the scattering coefficient: 



/3 



scattering intensity 

incident intensity-path length 



(207) 
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’.rhero tlie path length dx is that traversed by the incident radiation. 
Combining the above equations: 

=- 3 a MdxdA 

=> (5 = ( 203 ) 



Until nov7, r represented the radius of the radiation field and a the 
particle radius. Again, to be consistent vith the literature, r must be 
assigned to the particle radius. Using the new notation, if there exist 
a population of scattering centers each described by U(r), then the a- 

?rom all the ^articles of radius r between r and 



r +- dr is: 



(3 — cr ©, 0^^ c(A/ Cr) 



(209) 



The quantity dU(r) may be WTitten as: 



( 210 ) 



dU 



’.■7here ^ is the particle number density per particle radius. Then: 



<^0 = o- (^ ) dr 



( 211 ) 



(212) 



Let the lower and upper radii limits be r^ and V 2 respectively; vre 
have: , ,, 

In the case of o' then/3 - extinction coefficient 

for a volume of scattering particles defined by dV and equation (212) 
would represent the total energy loss by an incident beam passing 
through a volume of depth dx. Assuming an x (path length) dependence 
for , then the amount of intensity lost per unit incident inten- 

drex) 



= "/3(xjdx 

■X 



rex) 
XCx) = Ic exp 



(213) 

(214) 
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where I^is the incident intensity to the volume at x-0. Tliis last e- 
quation is loiovm as Bouguer's law. [23) It describes the degradation of 
the incident bean as it traverses a material path. In order to express 
the scattered intensity into a solid angle dii^we recall that the in- 
tensity scattered per particle into solid angle dH was given by equation 
(200), Let the parameter r in that equation be replaced by R as r now 
describes the particle radius. 

(215) 

As before, the number of particles in the cell, II, applied to equation 
(215) yields the total scattered intensity per imit radii per volume; 
then, for all particles of radius r between r^ and V 2 ' 

^ ^ I ( e) g C^) dr 

=> I = r ( nn, e, 0'^ 

For unpolarized radiation: 

c (ry), e) 3 (0) = ^ 

where i^ and i^ are given by equations (l65) - (I68) and (l70). For the 
polarized radiation: 

(-2 ( 218 ) 

Equations (212) and (2l6) are the specific results of the theory 
used to investigate aerosol parameters and v/ill be discussed in the next 
two sections of this paper. First, some comments on the collective 
scattering phenomenon are in order. It is v;ell knoum that the more 
ordered the scattering centers become, the more coherent the scattering 
due to the phase relations between sites. m This leads to a reduced 
baclcscatter u/ith corresponding increase in forward scattering. In 
crystalline substances, there is essentially no backscatter, when off- 



(216) 



(217) 
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resonance; and the incident beam merely undergoes propagation velocity 
changes. The condition met by a gas at atmospheric pressure or less, 
where the number density is low on a radiation v;avelength scale, ensures 
that the scattering is incoherent, allowing the summing of intensities 
from all the scattering centers. The Mie theory for a collection of 
particles requires a random distribution of scattering sites to achieve 
this incoherent intensity accumulation. Quantitatively, the near-field 
interactions between the particle may be ingored if the average distance 
between sites is at least ten times their radius. This criterion is 
satisfied even for very dense fogs and the independent scattering model 
computations are valid for practical tropospheric scattering. However, 
significant number densities ^ri_Ll lead to a large flux of previously 
scattered light in the region and add multiple scattering effects to the 
received signal, A quick checl: on the validity of the single scatterer 
model is to verify linearity between the observed effects, I and/3, and 
the number density of scattering particles. If the relation is non- 
linear, one must suspect the existence of coherent and/or multiple 
scattering. Alternatively, if the relationship is linear, it is certain 
that the independent, single scattering conditions are present, A 
measure of this condition may be obtained from equation (214), (2^ Let; 




(219) 



be the optical death of the medium, then: 

-r X 



X C^) = XqS 



( 220 ) 



The conditions are: 

rx<o.i (transmission > 9055 ) - single scattering predominates 
o.i^ Tx < 0.3 (transmission 75 to 95^) - multiple scattering 



effects may be removed via correction to the 
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single scatter model 



7'>i (transmission < 15 %) - the complete, non-ana- 
lytic, approximate multiple scattering theory 
applies 

Virtually all aerosols in the troposphere satisfy the first criterion. 

The literature abounds with material containing graphs of the Mie 
scattering functions for many specific cases. The Appendix outlines one 
such case for a specified set of conditions, Kowe'v^er, some general 
information can be reduced from this plethora of data. The amount of 
forv/ard scattering is about equal to that of the backs cattering at a 
size parameter (^) equal to 0,01, Asciis increased, the amount of for- 
v;ard scattering also increases with a liJce decrease in the backscatter, 

1, the forward scatter is on the order of 100 times that of the 
backscatter. Increasing still further, there e:d.sts a unique scattering 
maximum with the existence of smaller subsequent maxima and minima. 

The scattering function increases smoothly as (Rayleigh) up until 
/v/0,5. Thereafter, the structure is highly dependent on the complex 
index of refraction. The large oscillations in the behavior as®< is 
vaTied become smoothed as the absorption index increases. As the num- 
ber density, as a function of the particle radius, tai:9s on mors of a 
detailed structure, there is a further smoothing in the scattering 
functions as small differences tend to be averaged out. Conceptually, 
one has at hand either an actua]. particle distribution or a specified 
model of such. 

Given IT(r), for each r, we can compute the required I^('x',©) and sum 

j. 

the total intensities using 'l(r) as a weighting distribution. This 
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total l(e^, 0) then appears as the snoothed f’onction, 

Lilce the scattering function publications, there are also a large 
quantity of analytical nunber distribution nodeis in use. The verifi- 
cation and application of those nodeis has been the prinary purpose of 
□uch of the aerosol e:q:erinents conducted in the past 20 yeans, and the 
uor:: is continuing today. Is computers become faster, real tine inver- 
sion analysis ‘./ill lead to a more accurate determination of the "true'- 
number distributions, Tvo such models shall be mentioned here. Some 
natural aerosols may be described with a uuo-paraneter function utilising 
a multiplicative scaling coefficient, c, and a shaping parameter ,n), as 



nresenmed 



by Junge Q*] ; 







( 221 ) 



■ I I ~ ^ c! o o 



previously defined in (210) end may vary from a.bout 



2,2 to 4.0. As an example, for silicates in the O.l^m to 10.0 range, 
one may use; 

r, = O.OM/um , y V =v 3 

Another such model is in reality a set of models that includes the Junge 
model as a special case, Deirmendjian [l] has proposed: 



fCr) = 



^ -br 
~TZ — - Cir e. 



( 222 ) 



(note; the ^ used here is different from the size parameter.) The 
versatility of this approach is the presence of four adjustable parameters 
(a,®**, b, y) sJ-lov/ing for a vide range of experimental data to be fitted, 
Tl'iree specific cases of this model are; 



Haze C (continental): -P T 



- b'Tr" 



Haze II (maritime) : -p re 

o 4 -1.5 r 

Cloud : P^r'^e 



57 



Host of the other proposed nodal distributions utilise sone forra of 
equation (222) as their basis £und canned routines for their conputer 
modelling are readily available. 

T'.ro other effects that have been omitted in the present theory are 
Doppler shifts and Brillouin scattering. The presence of random 
(thermah) a.coustic vaves in an atmospheric cell will add a degree of 
order whereby Brillouin scattering of the incident 'wave can occur off 
of the resultant density fluctuations. However, this effect is extremely 
•.;ecL; for gaseous materials and the degree of coherence added to the 
structure of the gas is generally undetectable. Similarly, this motion 
of the gas will cause Doppler broadening of the wavelength of the in- 
cident radiation adding to the natural line ’width of the source, Tven 
for the sharpest sources, loners, the natural line'widths are on the 
order of several ilHs, ’.hercas, for almost all atmospheric conditions, 
the enoected Doppler shifts sne on the order of a few tenths of a Mlis. 
Therefore, the sliift is generally lost in the natural linevridth and the 
relative motion between the detector and the atmospheric cell may be 
neglected. 



III. LID.4R 

The sophistication of the laser systems available today, coupled 
’with the fast computer processing that can cheaply be obtained, has 
singled out the laser radar, LIDAR, system for increased use as a tropo- 
spheric e:<perimental device. The LIDAR should have the ability to per- 
form a 'wavelength scan with selectable band’width receivers. It is the 
atmospheric attenuation of the laser pulse and the introduction of 
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random phase and amplitude changes that combine to form the scattered 
pulse, or echo. The subsequent inversion of this data then allows the 
e:q3erimenter to assess the optical characteristics of the atmosphere 
on a real-time basis, ’.fith the system geared to perform azimuthal and 
elevation scanning as well as range gating, one can then map out the 
atmospheric optical behaviour within the visible hemisphere. The two 
most important atmospheric parameters to be considered for a pulsed 
LIDAR system are the extinction and scattering coefficients. The total 
extinction coefficient may be uritten; 



The other processes tahen together are typically two or more orders of 
magnitude weaJcer than the Rayleigh effect and tend to get lost in the 
signal noise, 

Rayleigh scattering is obtained from the Hie solution assuming «>c <<1 
using equations (184) and (1S9) - (191): 



Rayleigh scattering for air molecules assumes m is real and m^:^1; 



(<3 - -I- / 3 ^ ^ /3 A y- /^o 



(223) 



where R - Rayleigh (molecular) scattering 
M - Hie (aerosol + other) scattering 



A - Molecular and Aerosol absorption 
0 - Other processes (Raman, fluorescence, etc) 




(224) 







(225) 



where V = (4/3)fta^, the volume of one molecule. Let N - l/V, the 

s 

jUOXgCILL vinniTmcir* on o *f*.Vion • 




( 226 ) 
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where u = (6+3^)/(6-7(f) , the depolarization factor, required (227) 
to accoimt for the anisotropy of air nolecules, for which cT = 0,035. 
"he quantity 

afe d'i - Pi nun'll /"ij-p q-^-? 

(228) 



^(e,g^) = 0.0b07l (1 + 0.932 cos ^9) 



also accounts for the anisotropy. Let (°cu be the air nolecule number 
density (vice II), then: 



^ 229 ) 

where it is assumed all air molecules have approximately the same radius 
of 10 ^^m. 

The Mie extinction coefficient is dependent on particle radius 
through equation (212) whose only unknoi-m is dll(r)/dr. The unique sharp 
v;avelength feature of the LIDAR allows one to work well away from the 
characteristic absorption lines of the sample. Conversely, it also 
allows one to type the scatterer composition. For example; 

Laser 1: A.,°^ ^ absorption (resonance) line of a particular 



species; 



— > /v3| = /3p^i + (^/TU ^ (^jqj 



( 230 ) 



where is the resonant absorption extinction coefficient wiiich is 



strongly dependent on X ’with tabulated and 

= <^A 

_ f 

The is the number density of the species. 
Laser 2: ^^just off the resonance; 



(231) 






flZ 



Due to the relatively vealc X dependence for 
For an analysis cell a distance R away: 



( 232 ) 



(233) 



6o 



(2u) I(r) = 

and bade at the receiver: « 



Then; 



x; 



lo, e'^-^ 



•oZ 



Thus, the f®* for the particular species may, 
this equation using the inversion techniques 
For the typical magnitudes of the considered 



(234) 

(235) 

in principle, be found from 
outlined in the next section, 
cross sections, see Table II. 



TABLE II 



Process 


o' (cm / ster) 


Mie 


10-2’ 


10“^ 


Rayleigh 


10-27 




Raman 


0 

-P 

0 

1 

0 


10“^ 


Resonance Raman 


10“^^ 




Fluorescence 


4 10“’^ 





Similar to radar in operation, the LIDAR has a performance equation 
v:hereby the system characteristics and the appropriate atmospheric para- 
meters vhcih contribute to the signal are related. [23] Consider a 
pulsed LIDAR operating in monostaticrnqd.e (receiver/transmitter co-located) 
being used to analyze a specific atmospheric cell via its optical prop- 
erties. The problem geometry is outlined in figure 9. In this figure: 
hj. - cell height 
R - slajnt range to cell 

n -receiver solid aagle 
r 

- transmitter solid angle 

u 

L - cell length "gated" via receiver on-time 
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Monostatic laser radar 



/ 




laser 

transmitter 



laser 

receiver 



4 - 



( exaggerated) 



Figure 9 
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A - cell area seen by receiver 
r 

- cell area seen by transmitter 

The xollo'/ing derivation shall consist of counting photons at each 
stage along the beam path; converting to intensities has been previousl 
discussed. 

1 , I'lj. - number of ohotons transmitted oer oulse 

U X * 

Vj. - photon frequency (transmitted) 
hv>^ - energy’’ per ohoton 

U 

"xhV, - energy per pulse 

o "D - * 

- laser outout oouer 



- oulse width 



h Vt 



(236) 



2. for Goncstatic LIDAil, the scattering angle is IcO degrees (ff ) 
for backs ca,tter. Assume that all oackscattered ohotons are received: 



=> 






(237) 



The cell volume is determined by tlie Lim.'fl timing sequence and is: 

(233) 



The solid angle; 



n = A/R' 



V = LAta* 



(239; 

(240) 



The pulse width constrains the ma:-cimum system resolution by limiting the 
minimum cell length achievable: 

(241) 






where c is the speed of light in the medium. 

This condition assumes tha.t the pulse tail refloctixig off the front face 
cell does not interfere •ri.th the pulse front reflecting off the 
= 1 1 
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3 



. I;. - nunber oT photons incident on the cell 

(Sx - atnosnlieric extinction coefficient for the given 



transnittod 

N)u =W(e '-C%(s)d5 



(242) 



lio'/ever, to account for bean forning recuirenents, assume an optical 



system efficiency of 






(243) 



4. The baclcscattering coefficient nay be defined in accordance 
with equation (212): 

fo a(n)~ (3(1T) = , r,© = (244) 

where (’a’(Ti) is the usual notation. This equation represents the back- 
scattered intensity per unit incident intensity per unit solid angle per 
unit path length. The signal returned from the cell is due to this 
process alone. Bouguers’ law (214) states, for PaCn) assumed constant 



throughout cell length: 



X- e 



^ (ocr (J7)L 



(245) 



v/here I is that oorticn of the incident bean that nasses tlirough L un- 



affected. The scattered oortion is I = I - I : 

‘ s 0 






(246) 



where IT - the nunber of uhotons backscattered and a ra.tio of intensities 
s 

is the sane as the ratio of the nunber of photons. T^/pical tropospherie 
aerosols yield: 



fO crCn)L « I 



(247) 



(if - 0.02, overall error 1)'. ) 

The restriction this nlaces on L is an upper Unit of several hundrec 



neters; Thus: 

1.1 
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e ci (oa'(fr)L 



-> H' = A/l Per (ir)L 



(24G) 

(249) 



4'.e 'orinie on indicates the number of ohotons ocr solid anrie that have 



been oaclcscattered. Combining equations (243) end (249) 

^ Cf^)L e%p I'f ds J 



( 250 ) 



5. low tho echo must pass through the sane propagation path on its 
return to the receiver as it did on its nay to the cell. However, the 
extinction coefficient nay be different, (3^^, due to the accumulation of 
23rocesEes lil:e daman scattering, frillcuin scattering, Doppler sliift, 
etc., though individually negligible. 

- number of photons at receiver mirror 

(239), (250): 



(251) 



j (250): r r^r Tl , 1 

= Nt o-(n-) L fl . etc H L Wt M 






(252) 



6, The optical system vri.ll collect the echo photons received at an 
efficiency of ^ and convert them to an electrical signal via a detector 
v;ith a quantum efficiency o£^ : 



f = Nnn 



(253) 



where li^ - number of photons received and counted per transmitted pulse, 
due to bachscatter from the atmospheric cell of interest. 



Tnus : 



Mr = ^cr rtr)Lft.-^r'?q^e,%pi~I d^] 

I - I ■ ■ " ■ ' ^ ^ 



(254) 



This equation nay be recast in the form of intensities if the photon 
population returned is large enough that individual photons need not be 
counted. Let D be the optical system aperaturc (D, - transmitter; D - 

"0 X* 

receiver), then the intensity becomes: 
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(255) 



I = NK v>/d 

Re^vfriting equation (254): 

Ir 



•iting equation (254): r r'Vu./x ^ , C] 1 

■It fe)(37) 7* ^ 



( 256 ) 



Thus, the important parameters for the LIDAR analysis of an atmospheric 
cell are the scatterer number densities and associated cross sections. 

7or a given system design, one has available the quantities it>%> 
A^, D^, D^, and can select the values for , L, and R, If one 
assumes a slow moving target cell, then the pulse round trip time, t^, is 



tp = 2R/c 

R.= Cc./0^r 



(257) 

( 258 ) 



Thus, one determines the range, R, by turning the receiver optics on a 
time, t , after the oulse is emitted. In order to avoid the oroblem of 
range ambiguity, lilce radar, one must ensure that t^ be less than the 

X 

period between pulses: 

tr < l/PRF ( 250 ) 

v/here PR? - pulse repetition frequency 
The unambiguous range limit may be written; 

= (c/e)/PRF (260) 

v;here R^ - unambiguous range 

The cell length, L, is determined by the time that the receiver is turned 
off again a time At after t^. For the entire cell to be scanned, the 
pulse must traverse L tv/ice in At; 

L= Cc/2)At (261) 

and the receiver must be turned off at time; 

t - tp t A t (262) 

One can then select L and R by programming a timing sequence for At and 
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t^. Finally, the interpulse period, T, which is a system deadtime be- 
tween pulses available for processing, becomes: 

T. (I/PRF) - r* (263) 

Now one measures and ~))f and applies equation (256) to determine the 
desired quantities ^’it) and This inversion is covered in detail 
in the next section. Note that when the equation for CJ from 
(202) is used to replace the in (244), then equation (256) 

for has the correct form as given by ( 21 6 ) for the scattered 
intensities. 

Thus, a model has been constructed ’whereby the system parameters are 
related to the cell unlmowns in a definite manner. But first, some con- 
sideration must be given to system noise. During the time the receiver 
gate is open, ^ t, the system vri.ll collect and count the background sky 
photons as ijell as the signal. The continuous ’’dark current” in the 
detector yields an electrical noise as a result of dark counts. These 
counts arise from photo electro ns being emitted from the detector photo- 
cathode surface even v:ithout an illumination present. Each of these 
noise sources may be reduced by proper design of the equipment and ex- 
periment. In the case of a detector, its dark current, I^, is either 
specified or may be measured in the laboratory: 

^ ( 264 ) 

v>?here - number of dark counts per sec 

e - electron charge per coimt 
G - detector count gain 

The number of dark counts detected as part of the return signal in a 
time At is and : 



/ r-; 

O / 



Mj) - (265) 

Let N„ represent the background count rate such that lUAt is the nun- 

D " d 

ber of background counts received per laser pulse. Assume a baclcground 
spectral brightness, that accounts for all the natural and artificial 
baclcground sources. This V/-^ represents a power density per solid angle 
in the background sky. Then: 

A \ (266) 

h V 

where ~ previously defined 

V - frequency of background photons (at A) 

A A - receiver bandwith, centered about 
Let II^ be the total number of noise counts per laser pulse; using 



equation (26l) for At: 

Mr, = (2L/c)(Mj,tNs) (267) 

Due to the pulse repitition frequency, PRF, the noise count rate: 

(PRF)Nn 

v/here N - noise count rate 

The statistical nature of noise yields an effective RIIS noise count rate 

of '/if. The signal has a count rate of: 

Sr = (PRF)Nr (269) 

v;here S - count rate from signal echo alone 
r 

Thus, the total counting rate is: 

S = PRF (Nr N) (270) 

Talcing the RII3 of the signal plus noise, one can form the signal-to- 
noise ratio: 

siH = 5r/^^3- (271) 

Assuming a one second time averaging: 
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S/N = V PRF 



(272) 



If the noise average, h , is not accurately neasured, the '.I. is replaced 

n n 

by 211^ (measured) in the above equa.tion. hote that: 

1. li! 0 => (s/A/)o ^ VCPRF;Np (273) 

2. if Nn= - C^I^)(i/\rr (274) 

Therefore, the S/h ~ PR? and by increasing PR?, one increases the signal- 
to-noise ratio, Kov;ever, via equation (260), the unambiguous range de- 
creases faster than the increase in S/P. Since a lange S/li is desired 
•;ith good raaige resolution, one selects a detector v/ith a lou j^ 7 ,(-'Q) and 
designs the system for lov; background reception (P^). If the time aver- 
aging had been porformed over a period of T seconds such that PR? T 



pulses had been transmitted and received: 

S/ri = VPRFT i^ir/V^/r'^A/n 

Ponbining equations (265)- (267) and (275): 



N = 



sTpRF r A/r 






(275) 



(276) 



Inamination of the above equa.tion and ( 254 ) leads one to some general 
system noise qucP-ifications. 

1. One desires a lo'.f detector I^ and a high'^^. 

2. Select a nonrow receiver spectral banduldth, A X ; another reason 
for the choice of a laser. 

3. Reducing L reduces the noiso as uell, but the transmitter must 
put out greater pouci- to maintain the sane signal. Then, for the sane 
average po'ver, the PR? must be lowered and the effect in S/P is only 

\| PPIT as cempaned to the 1 decrease and Pj. increase, increasing S/'", 

In general, high po’-.'er/lov; PPl yields a better S/P than a lov/ pov/er/ 
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4. A snail requires a louer bean divergence at the transnitter 
and on optical matching 'oetx-ieen the receiver and the echo bean pattern 
(another laser capabilit;/) . 

5. Without any noise, there still exists a finite (S/:I)q vliich 
represents the "shot noise" statistical nature of the signal. It must be 
enphasised that all the quantities in equation (276) represent the nesn 



os tne .inanccvoec: s'craaoj-os. 



6 . Mr « Mr =<> S/N ~ Nr 

If Mr Nn => S/^J ^ V777 

The actual tropospheric neasurenent will include both 
(noleculcn) and hie (aerosol) scattering. As an exanple: 
a. single ’.rater droplet, r = 1 (hie) 



:277) 



(2-76) 



lO”® cm 
(^ '^ 1 0 cm 



-3 



(3 



m 



10 



-7 -I 



cm 






"" 10 



*2& E 

cm 



=> 0R 



rs -3 

10 Cm 






Thus, the addition (3j^ to the total (3^ is on the order of a far; per- 
cent and any quantitative neanurenent must t?J;e into a.ccount both the 
aerosol and noleculan scattering. The t’virbidity is defined as the ratio 
of the total scattering intensity to the Rayleigh scattering intensity: 

= It/Tr 

oone researchers attempt the neasurenent of this quaaitity independent of 
tile intensity, then use the res’ult to deduce the aerosol I... .Another 
nethod assumes that the a.erosol concentration above a certain height 
(e.g. ^ 30 I:m) is sufficiently lov; to allou it to be noglectcd. Then 
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( 230 ) 



equation (256): 

' Ir ^ 



[Text, one computes from (256) assuming the conditions of the U. S. 

Standard Atmosphere and equation (226). ThuLs is the e:cpected echo 

intensity due to molecular scattering alone. The ratio betvreen the mea- 

sured and the expected yields a scale factor, b, that corrects the 

standard atmosphere for the local conditions. Using this factor, one 

assumes a density distribution for the local atmosphere the same as that 

for the standard and commutes the e:cnected value of i_ at the cell's 

R 

location. This intensity is then scaled using b. The measured value of 
I_ is then corrected for this Ip(b) and the result becomes I,., the Mie 



intensity: 



(2S1) 



For other means of separating the molecular and aerosol intensities from 
a measurement, the reader is referred to the v;ork by Hall, et. al, Lvl 
This Rayleigh-atmosphere assumption is the starting point for the in- 
version process. 

The Standard Atmosphere may be represented via the folloving equations 

temperature T('*F) = 59.31 - 0.00357 (2S2) 

pressure P(a4m)= i. 022. 157/ (2S3) 

Mass density _ O,002M25e/p (-3.3^5^ (284) 

(slugs/ft-^) 

where h may vary fix)m 0 to 37,000 feet.[)3j The more precise form for 
these relations are also represented as canned computer routines. 

The critical term in equation (256) is the where the sub- 

script M refers to the aerosol distribution alone. The extinction com- 
ponent in the exponent of e in this equation represents the scattering 
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and absorption cross sections integrated over the entire 4'ff steradian 
sphere. If one assumes some representative model for the aerosol dis- 
tribution, as was done for the Rayleigh case, then the only unlaiovm in 
(256) is The two averaging processes, integrating aver 4 'n' and 

integrating along the path R, reduces the sensitivity of the exponential 
term in (256) on the actual aerosol distribution. After completing the 
inversion process, one may input the best estimate for tliis distribution 
into the exponential term in (256) and then repeat the inversion, Ko’j- 
ever, the processing time involved may prove to be the final constraint. 
The LIDAR experiment yields the scattered Hie intensity from the 
aerosol population in the analyzed atmospheric cell as a function of 
the pre-selected parpmeters at the laser (e,g, wavelength and polarization 
scanning) , '.7ith this data as input, the system is prepared to perform 
the inversion routine. Some computational needs are reduced for the case 
of laser backscatter as the scattering angle remains fixed at ISO degrees. 
This baclcscattering condition simplifies several calculations as out- 
lined in the Appendix, 



IV, MATHEMATICAL PR/ERSIOIT 

Tlie assumption made throughout this paper implies that the e:<peri- 
menter need only apply the exact Mie theory to his observations and per- 
form a simple data inversion to discern the aerosol properties of 
interest. However, this inversion is by no means simple. Aerosol prop- 
erties may be separated into two distinct catagories. The first con- 
siders those parameters that are explicit in the previously defined 
scattering functions; i,e,, spherical aerosols, aerosol sizes, low 
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number densities, and the relative indices of refraction. Then there 
are those aerosol characteristics that require some approximations (in 
many cases, severe ones) to the theory; i.e, , non-spherical shapes, 
anisotropy of the shapes with respect to the scattering plane, and high 
number densities leading to mioltiple scattering effects. The first few 
LIDM experiments investigated the actual aerosol ranges for these para- 
meters, They found that in excess of 90% of the real aerosols studied, 
the spherically symmetric single scattering theory was valid under favor- 
able atmospheric conditions. These conditions require the absence of 
dense water clouds in the atmospheric cell and a depolarization of the 
incident beam limited to a ma:d.mum of 2%, with 1.5 (j of less actually 
measured. [^ 3 ] In the following analysis, these favorable conditions 
will be assumed, as will the requirements of single particle scattering. 
These conditions are enhanced by the common field of view shared by the 
LIDAR transmitter and receiver. 

Recall equation (21 6) for the scattered intensity a distance R av/ay 
from a cell containing an aerosol population described by di/dr with 
radii limits r, to r^: - / ui, 

,=.« 1 144. " 

The first simplification is to replace the particle radii continuum with 
a set of distinct particle classes each with a characteristic radius, 
number density and complex index of refraction. Thus: 

r = X(rripm2.) n > ''iv ') (235) 



where m. - j 



th 



oarticle refractive index 



.th 



r^. - j particle radius 

N. - particle number density 
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X - incident radiation v/avelength. 

^ - scattering angles 

The large number of unknovnis that need to be determined will ensure a 
non- unique inverted solution from the limited sample of observations-.a- 
vailable. There is also the added difficulty of variance in the sup- 
posedly Icnovni data v/hose uncertainties at best generate unbounded noise 
in the final results. Therefore, some attempt must be made at reducing 
the initially large number of unlmno’.,Tis. 

The dependence of the intensity on the size parameter, has been 
discussed previously. Briefly, the non-monotonic fluctuations as is 
varied is the distinctive feature of the Mie theory that enables the 
process to be of good use in tropospheric ar.alysis. These ma::ima and 
minima are reduced by increasing the field of view and/or the transmitte 
bandwidth. It is to LID^iR' s advantage, then, to retain this complcz 
structure as a data filter. See the Appendix for a sanple plot of 

The complex index of refraction, n, has a dispersion relation lilce 
any electromagnetic medium. See Table III for an example. 

The less than 1/ change in the index of refraction is t’-grical for 
dielectric media. The change i:' the absorption, y" , is even less, otudi 
have indicated that ti.e dispersion in m is balanced by the randomness of 
the aerosol size distribution in such a manner as to almost negate this 
dispersion, 

The aaigialar dependence is a pre-specified set of conditions, espec- 
ially for the case of a LID^IP. oachscatter system. This system consists 
of a linearly polarised '.rave undergoing only bachscatter observations. 

As a result of these and similar considerations, one can reduce the 
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4 









number of unl-uiovr.s drastically by malcing onj.y a feu simplifying 
assumptions. These ha,ve also been verified by the laboratory ejialyses 
of many aerosol samples, natural and artificial. 

TA3L1 III 

Tjicuid water (Q20°C) index dispersion; m = n('!-i')C) 

(A°) n 



4046 . 6 

4471.5 

501 5.7 

5460.7 

5292.6 

6562.3 

7065.2 



1.342724 

1.339423 
1 .336363 
1.334466 

1.332963 

1.331151 

1.33CC19 



1. The measured values for n raiige from 1.33 to 1.5 



VI 7, n 



y .3. 



G 7 



to 7.18 surrounciin,2 oar. 
be cboro.ctGrised by: 



lass cf aerosols (a oorticulo^tcs) ray 



no = I. 5 4 ~ lO. Ooi) 



' 2... O y 



Thus, assume that the complex index of refraction is both hnox/n and 
unique. Only the pas’ticle’s radius arid number density is alloi/ed to vary. 

2, T;rpical non-precipitating aerosols in the troposphere are described 
with their radii between O.Olyxn and 15ykm. Their number densities can 
be represented by distributions containing only a fe’.; parameters. The 
important model distributions \/ero presented eanlier and the usual, goal 
of an invers3.on is to best fit one of them. The initial estimates, 

\;here required, of paarticulate properties usually includes such models. 

The set of unlmoams reduces to the particle number densities ci'xl 
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their radii. 3y the use of pre-selected values for r, one only has to 



solve for the distribution functions granting that the radii range is 
properly considered. Therefore the paraaeter of prinary interest is the 
scatterer nuaber density. 

Define the ilie function, K: 

(2S7) 



Ulit: J. Uliu U 



and the particle size distribution function, f(r): 






dhiCrl 
dr 



lie scattering intensity, equation (21 6), becomes: 

I M = r V f 



(2SS) 



(2S9) 



where x represents the set of user defined parameters (m, X,0, 0). The 
value l(x) is measured, the quantity K(x,r) is computed and the unlaioun 
f(r) must be solved for, Tliis relation is classified as a Tredholm in- 
tegral equation of the first kind. The kernel of the equation is K(o:, r) 
and I(x) represents the transform of f(r) utilizing the stated kernel. 

There exists only a finite set of measurements, l(x), over which the 

'th. 

inversion may be performed. Let the subscript i represent the i 
measurement of a set of n observations: 
luo = = = / K (>^i,r) f(r) dr (290) 

' I 

The notation, g^, has been introduced and v/ill be used later in the paper. 
The inversion process is by nature an unstable one and this may easily be 



demonstrated . 



Define: 



sinFr] dr 
Jo 



(291) 



where i = 1 , 2, 



n 



F- 1, 2, 3, ... 



and C is an arbitrary constant. 

(290): 



jetting r^ = 0 and ^2 - in equation 
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r It 

‘ ^ c sin (Fr) K (^i ^r)cir (292) 

Since the sin(?r) acts as a weighting function for the kernel on the 
interval (0,11), as ? increases '.-dthout bound, the second integral tends 
to 0. Thus: 

One can anoroxinate the transfomed values of f v;ith the function S_, as 
close as possible by talcing F large enough. Yet the josinYrj function 
nay be nade arbitrarily large and 

?(f) = p(r) + CsLnFr (294) 

will transfom as f(r) with little to no resenblance to f(r) for large 
values of 0, Therefore, one can generate transfoms that cone as close 
to the desired values, g(x), vdthin any desired accuracy as one may vdsh 
and still have no idea v;hat the accuracy with respect to f(r) nay be. 

The technique used to remove this instability is to require the ad- 
dition of a constraining relation on f(r) other than equation (290). 

Tlois nay be done either e:q)licitly or inplicitly. The proper relation 
must contain the desired qualities of f(r) requiring a high degree of 
a priori Icnouledge. Out of the near- infinite set of possible solutions 
for f(r) generated by the inversion, the constraining relation v;ill act 
as a filter to obtain the unique, or near-unique, solution. 

The presence of measurement error in the g(x^), larger than some 
quadrature error, allows the use of an appropriate quadrature approx- 
imation of the form: ^ 

>3 ^ ^ 

for k intervals of the argument y and the w. are the quadrature 

o 

weighting coefficients with increments Ay.; . For equation (290): 
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( 296 ) 



X f (rO A r] 

0-1 ** 

v/here i represents a aeasurenent and j a radii interval. 

Define; for n measurements: 

1, column vector G: i - •••> (297) 

2, coefficient matrix A: a. . = v. .K(x.,r.) (29S) 

ij ij ' 

j = 1,.,., Ic; Ic radii intervals 

3, column vector II: II. = f(r.) Ar. (299) 

J JO 

Then, equation (296) may he :n:itten as a matrix multiplication: 

G = .III (300) 

The value N. represents the total number of oarticles per unit volume 
J 

’.'.dth radii between r. and r. + ^r., as can be seen via (233): 

- J^. f(r)3r = fCrj) [( D + A rj) - rj] (301) 

and f(r.) is considered constant over the inter/al J. It is further 
J 

assumed that the scattering functions have only small variations with- 
in each sub-inter'/al radii range. 

The problem has now been reduced to solving for the k unlcno’.ms in 
II given the (k)x(n) computed values in A and the n measurements in G, 
It must be emphasised that a quadrature approximation v;ill introduce 
serious error in f(r) unless the quadrature v/eighting formula is of an 
accuracy higher than that of the expected measurement error. Let Hq 
be the nu,mber of independent observations made such that n., ^ n. If 
k = Uq, then II may be solved for exactly: 

(require A to be non- singular) 

=r> A''A = i (302) 

I in a matrix equation represents the identity matrix, 

A^J = G => M - A“'G , (k-Ao') 



where: 

( 300 ): 



(303) 



A 

where N represents the interpolated solution vector. However, if k 
n^, then M may be solved by minimizing the square of the error (least- 
squares) : 

(require a'*A be non-singular s> A*' ; transpose of A) 

(300): (304) 

this step forms a square matrix which may be inverted, 

M ) ' A y ( K n J (305) 



That this last equation actually minimises the error may be visualised 
by considering an appropriate spatial geometry which contains the vector 
G and the vector AH^ where N. See figure 10. The vector C— is 

seen to be representative of the error vector betv;een G and AH^. The 
minimum error occurs when this error vector has the smallest magnitude, 
thus : 

(AN,) -(<--flN,) =0 (306, 

(the • represents the vector dot product). 

A 

Then for any AH in the range of A, the vector A2I^ is still perpendicular 
to it; 

(^^j)*(G~AAi) =o (307) 



Writing the dot product as a matrix multiolication requires; 

(flW,)^(G- a 4) = o 

And uy r 0 since it represents any vector in the range, thus; 

(304) 



( 308 ) 

(309) 



In the work that follows, when vnriting the inverse matrix A \ if the 
matrix is not square, the quantity (A.'*’a) ^A"*" is implied. 

The procedure as indicated by equations (303) and (305) are in large 
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n-3pace vector geometry 



G 




Figure 1 0 



SO 



oone 



unsuccessful as they do not contain the constraining criterion. 



further problems concern the size of k. 

A 

1, Small k: the vector II does not accurately represent II due 

to the r. set sparsity and/or the crudeness of the quadrature. 

J 

A 

2. Large k; the components of II taice on an oscillatory be- 
haviour irrespective of IT and the coefficients in A become particularly 
sensitive to round off error. 

This sensitivity in coefficients is best characterized by the con- 
dition number of A (or any matrin) [5] : 

conc/(A) = jflj lf\''l (310) 

■!hers )A| - xiora of a matrix 

The norm of a matrix may be defined in any one of three uays depending 
upon the desire of the user to highlight specific as'oects: 

(311) 



1 fli = 


(Y10.X 


1 . . 1 




1 ' ^ 1 1 


1 $ J ^ K *(. > 1 


kcjl 


: niaxirrtum ctLum\j - sum 



v/nere 



a. - magnituae oi tne component a. . 
^ 3 ' ° K. * 



: maz-timum ro-..-sum (312) 



% 



I I CO = I i n I ^ ^ 

(normaj-ly used) | /1 1 g =• ( jn | ^ ) *" ‘ Lnclidean norm (313) 

Lecause A~ ' must be computed and '«ill contain at least the inaccuracies 
in A, the cond(A) vdll be costly to compute and inaccurate as veil. 
Ilov/ever, the test is th.c magnitude of the cond(A). If large, the matrix 
A is termed poorly- conditioned and the lil-celihood of a successful inver- 
sion decreases rapidly to zero. Corrections, lilie scaling or partial 
pivoting, must then be applied to A. Ilote, for the identity matrix, 
cond( 1 ) s 1 ; 

cond (a) > I (314) 

since must exist. Other oromerties of the norm are: 



Si 



lfl[ ^ o 

lf\i - O /? = O 
1/3/1 I = /3 //? I , K? > O Oind a. nurvher- 
l/l+Bj ^ Ifll +IbI 
1 B/1 I ^ 161 lAl 



(315) 

( 316 ) 
(317) 
(31S) 
(319) 



The ccnd(A) is nost useful in bounding the error expected in the computed 
solution and relating this error to the nagnitude of the residuals. 
Define: 

.hi = 3 (320) 

’.'.'here X - exact solution vector 

A 

X - approximate solution *'ector 
3 - vector of xno’..Ti quantities 
A - operator matrix 
The residuals are: 

where R - residual vector 

A 

^ r: - T 



A 

Ax 



F> = 



(321) 

( 322 ) 
(323) 



using 3 bp inverting equation (323) end the norm relations (31 5)- (319), 
one may incorpora.te the cond(A) as an error measure: 



4 



Rl 



CoMd(A) 161 IXl (324) 

A 

Therefore, the relative error in X, Gonta,ined in 3, can be as great as 
the relative rosiduaJ. times the cond(A) and as SDcill as this sane re- 
sidual divided by cond(.t). For large cond(A), the residual p-ields little 

A 

to no information on the accuracy in X. Xquations (320)-(323) may be 
used to define an iterative improvement algorithm (IIA). 
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« 



r 











1 



A A 

. compute X by = A 



(325; 



(321 ) 2. compute R by R - E - i-Z 



3. compute 2 by 



_ ,-1. 



( 326 ) 

^ A A 

4. correct X by X = R + X (327) 

If the system matrix A has errors in its coefficients, the system 



may then become ill-conditioned even if the residuals are small. 



A 

A = A + ?: 



(32S) 



i.'here A - true coefficients 

H - error in coefficients 



The equation actually being solved is: 

AA 

A ’ -r "D 



(329) 



(330) 



'7nen, analogous to (324) : 

iifrSf i cowt> (a) IHt 

Ix' 1A\ 

The relative error in the residual can then be no larger than the rela- 
tive error in the coefficients in A times the cond(A). 

Finally, an approximate accuracy check may be obtained from [p3 : 

1 E or |-( I _ p , , 

j -r (33^) 









-p 

lO , 



then X is probably correct to p digits. 

Equation (300) is assumed to be accurate within a small measurement 
error and equations (303) and (305) are considered insufficient to 
remove the inversion instability. There are several techniques avail- 
able to remove thses instabilities of which tliree will be briefly 
covered. 

The first method attempts to smooth the instabilities by one of two 
means. Define: 

u(w)=((r-AM)\6-ANi) ^ y, (332) 
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where u(il) - the smoothing measure (a number) 

^1 ’ ^2 ^ 0, numbers 

Mq - a priori information function, vector 
B - matrix describing some desirable smoothing of N 
(G - Klf (G-AII) - a measure of the accuracy with which N satis- 
fies equation (300) 

"ip 

(ri - IIq) (N - ITq) - a measure of the departure of N from the 
expected 11^ 

(bn)''* (bn) - a value that measures the departure of N from ideal 
smoothness (BN = O) 

The smoothing technique involves the minimizing of u(N). Some matrix 
differentiation rules will be needed 0-3 : 





= zx 




(333) 


AAX 


= A 










(334) 






AZ. 


(335) 


<iX^ 







where X, Y, Z, are vectors and A is a matrix, all of the proper dimensions 
such that the indicated products are defined. 

Condition 1 : ^2 - 0; requires a priori knov/ledge, Nq. 



(332): u.Cw)= 




(336) 


where VoC^) " 




(337) 






(338) 


The required condition is: 






= mimraum => U.C’^) 


= 0. 
N 


(339) 


Utilizing (333)'-('335) ; where 1 is the identity matrix: 




^ = - 2/^’'Cc--an) 




( 340 ) 






(341) 



-4 



Use (340) and (341 ) in (339): 

N, = 0 ■'■^1 ) (342) 

This is the required solution to equation (300) given the constraint 
function Nq. Note that, by setting = 0, one obtains equation (305). 
Since the multiplier, Vi , represents the difference magnitude bet;;een N^ 
and Nq, one might e;<pect to use the estimated error, (331), as a first 
guess at y, . 

Condition 2: =0; requires no pre-conceived solution. 

(332): IIjCn) = y, (W) 3/z YzM (343) 

v;here y,, - as before (337) 

(344) 

The most common selection for B requires the solution to be smooth 
through the second derivative. In difference equation form: 



^ dNj ^ ^ Ki' 

Nj = ar," Arj = rjn -r7 ' 

^ A Nj ^ N'j*i J /Nju ~ * 1 Mj+i ” Alj 

'i ^ dr -4 rj ~ rjti-rj ~Arj Ar j..i A rj 

Assume equally spaced quadrature points: 

^ CK-I ) 

Tj r ^ Ar 



K - I 



Note thatAr. is independent of j. Thus: 

< =iv*(Mi- 2 A/iH * Mj,x) 



(345) 

(346) 

(347) 

(348) 

(349) 

(350) 



Since only the coefficients of N. are important, ignore the l/Ar scale 

J 

factor and vo*ite: 



N" = ( I 



-z 



0 



A 

Mj*i 



85 



(351) 



or II - 31! 

in matrix form where 3 is a (l:)x(k) matrix defined as; 

/, 



(352) 



1-210 



0 

0 



1 -2 1 
0 1 -2 



B = 



0 0 0 1 



0 0 0 
0 0 0 
0 0 0 
0 0 0 






0 0 0 0 
0 0 0 0 
0 0 0 0 



1 -2 1 

0 1 -2 
0 0 1 



(353) 



J 



The vector 3II is thus the discrete analog of the second derivative of the 
unlcnoun vector II. Other forms for 3 may be used, notably reference (^2(fl . 
As before, the constraint shall be to minimise the quantity U2(H) In an 
effort to achieve the smoothest function II v;ith a minimum of error in GII: 



U 2 M- minimum => ^ j ^ =0 

and Y 2 : 

How, (340) and (355) in (354): 



(354) 

(355) 

(356) 



This is the required solution to (300) v/hich satisfies the constraints 
outlined above. Again, setting ” 6 yields (305). A recent work 
into this particular solution has yielded considerations for the choice 
oi if the factor ^('2 is too small, the instabilities arc not 

removed from the final solution. And if chosen too large, the system 
becomes over constrained and indenendent of the measurements. The best 
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choice appears tc be ^ sr.d should be so chosen until further 

efforts provide a better choice. 

Both solutions via this technique and the results from the next 
method nay be improved in the iterative fashion as outlined in equations 
(321), (325)-(327); IIA. The algorithm: 

(initial estimate) 1, obtain by (342), (356), or (366) 

2. - G - (357) 

3. (35S) 

4. = v(°) + y(‘'^) (359) 

5. is < some small positive number ^ ? 



(m^^‘ residual) 

(m * error estimate) 
(nev; estimate) 

(test) 



no: m = m + 1, go to step 2. 

A(rn+1 ) 

yes: done, the solution is h to an accuracy on the order 



of 5*, 



The next method involves the computation of some statistical para- 
meters. Only the general outline shall be described. See reference [2l] 
for details. Let and be positive definite matrices (all eigen- 
values positive) that act as veighting functions, then: 



v(f^) = (G-AM)^ri + C/^-/Vo) 

uhere /7 is (n)x(n) and /Q is (k)x(I:). 



( 360 ) 



Let IJq be the mean vector for II as compiled from previous statistical 
data (model). Define the expected value operator E as: 

l,rna g (X-,) -Pro ( 351 ) 

where f(x) is the probability density function for g(x). Let the vector 
£ have zero mean and measure the error in G. Define the covariance 
matrices thus: 
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(363) 



5 . = E&E-3 



Assuning a nornal distribution for 6 and 'I, then the Gauss-Maricov 

A 

theoren yields a statistically optiaun estimate II by minimizing v(Il) 

s 



for: 



r, - 5, 



-I 



(364) 

(365) 

( 366 ) 

(367) 

A _1 

The resulting covariance matri:: for (II - II ) is found to be X . Thus. 

s 



then : 
where 



G = 

X = sI'*a"SI'a 



the expected mean square error may be written as: 

= trace X 



(363) 



Without solving equation (300), one can thus obtain an accuracy estimate 

A 

for II under the stated statistical conditions. One can also annly the 

IIA, (357) - (359), to enhance the estimate II . 

s 

The last teclinique to be considered is an iterative scheme somewhat 
analogous to the iterative Improvement algorithm. The method used is 



the Landweber iteration . [l . Ir. general, an iterative guess 



.( 0 ) 



made ba,sed on some ?, priori Icnowledge and then improved successively 
to form a sequence w'hose member, I. ‘‘‘ , transforms to the measured 



( m ] 

' for smaj.l m. Tills 



G as m tends to infinite. The techniques' success is due to the auto- 
ma.tic constraining of the oscillatory nature of II* 
ellcwo an approximation to C- early enough to avoid the oscillations if 
the secuence is terminated when the residual error is on the order of th 



error in G, TIio Landvreber algorithm uses; 



(369) 



where reflects the eic'ected solution. If the coefficient mc.tri;: 



S3 



l:as "j-l _?ositiv9 con_,onent3, the conver^^ence i/ill 
finins e. diac-ozial :3atri:c D: 



D: dji - ( ^ ' 



I** I 



J = I,..., K 



eXid reva’itir." the estimate as: 



Cm) A Cl"'’) «T/'/- nu 

=N^ + tiA’^CG-flWo ) 



be ir."ro'/ed b’^ ds- 






(371) 



This conpletes this short survey of inversion techniques. A coanon 
feature for ell of then is their use of a priori loiouledge to constrain 
the selection of probable solutions in order to enhance the choice of 
the proper solution. Indeed, this criterion nust be net by any tech- 
nique the ezroerinenter chooses for inversion, if he is to have any 
confidence in the results. 



This paper has endeavored to present the requirenents for a device 
designed to nap the tropospheric aerosol properties on a real-tino basi 
The specific ncdel equations ’..-ere derived alter a theoretical devclopne 
\ic.3 perforned in such detail tha.t the user becane faniliar ’.:ith bhe 
linitations of the theory. A candidate for such a device that uouj.d 
utilise advantageously the specific theory characteristics is the la.ser 
radca-, or lIDfl’,. The LID.'l syston allo\<s for the direct selection of 



nuts for the data inversion to bevin. A feu Inversion technloues were 



subject was by iio neans exhausted. The aerosol LIDri?. nust be capable 
of polanisation neasurenents and fast wavelength scenming to increase 
the nunber of independent observations nade. 3uch systen canability 



S9 



aecrease 



\/ould 



:3ot 



_asers 



aerosol 



th.e rcouired enour.t of 
r.Uot decide on. 

h nresont ond future developnen 
argue strongly for the design, 
analysis system. 



a priori laior.'iedge the e::peri- 



ts in fast processors and versat 
constriaction, and use of a LIDif? 
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ID SiQCRITid'I 




This section presents 


some 


1 further equations to 


facilitate the 


compu- 


tation 


of the required fuc 


itions from the Mie theory. 




Bachs c; 


itter: 0 = 


lc0° 












rr„ (i 80 ’j 


— 






(372) 




TnOsoh 


- 


'TT, 


CP 

o 




(373) 




5, (l90'} 


— 


1 


C-O'^Cn + x )(^n -An) 




(374) 




S^(l 80 ’) 




-s, 


CP 

o 




(375) 


bac 


:kscatter g; 


lin: 


6 




1 (-0 Cn-*2)(b^-o-n') 


1^(376) 



Scattering functions : 



Scattering functions : 

recurrence relations : 

7^ (cos e) = C 05 © Tin ©) - sin^eiTn (cos e) (377) 

ir„ (C 05 e) = (cm eX^) ®) (373) 

■n'n(;6ti5o)=f2'’-')rr„.,(“Se)+ft'l2(co5 e) (379) 

1("o (C05 ©) = o (3S0) 

TT, (cos 0) = I (311) 

i?;'(wse) = o (302) 

1T'/(coS©)= o (303) 



Scattering coefficients ; The equations for and may be vn^itten in 
terms of the logarithmic decrement function, (Z) . The equations 
below have been so x-rritten and then incorporated into an algorithm that 
was used to compute the scattering and extinction efficiencies for a 
particular case of m = 1 ,54(l-i0«C0l ), These values vfere then plotted 
on the enclosed graphs to picture the oscillatory behaviour of the 
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scattering efficiency, Q , for <=^ fron 0»1 to 17.0 in 0.1 increments. 
The data points were connected by straight line segments to complete the 
visualisation. It must be realised that further structure may appear 
under finer resolution. The routine was programmed for the heulett- 
Pachard model 410 machine vrith the program listing accompanying the 



figures. 





(384) 


m : ^eptC'- ■<•><') = 


(385) 


/3 - mcx = /0-< 


(3S6) 


llote in the follouiig that orJLy ^ requires 


a complex argument. 


(iSo) QetiM - -ij |_^(4n+0 

(187) QscaC'*^) - iy + 1 1 






(387) 


, m 1 c/3 )_ 1 

'='> -■ LT“(»0- -”7? J 


(383) 


9i C°^) = sin 


(389) 




(390) 




(391) 


Xo C°<) = C05 C^) 


(392) 


X, M + % (■.=<) 


(393) 




(394) 


(1U) tW+iXn(°^'' 






(395) 


x;;m= 


(396) 




(397) 


i'V) - 3n 


(398) 

(399) 
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/ 



(400) 



tar? /^* 



.Al,?:crith.'a 



(2n.0[|a„(=<)(4|t„r-‘)rj 
=> Qsca(f“) = ^ %' (^) 

<Jz(«)= (2n ♦ i) f?c fa, C«<) H- t„ C»<) j 

=> QextW z: ■^% 2 .(^') 

Begin: A, 1, input: opt -i ^ 

2. coapute: rrtx ^ y 

B. 1. input: (S^cT (required accuracy) 

2. compute: 5 y 

3. assign: /7r/ = ^2 = O 

4. compute: y, , y, C^) 

KoC^^) 5 X, C»<) 

(0) 

C. 1. compute: ("oc ) ^ ^6>< ) 

C”-) 



2, compute: 

3. compute: 



7», r'i^) 

(^) ^n(^) 



4. test: is 



> <§ ? 
C*=<) 

yes: go to D.1, 

no: go to C.5. 

5. test: is ^ ? 



^ C“) 
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yes: go to D.l. 



no; go to S.1. 



D. 1. 


assigTx: 










2. 


compute : 


(^), C°<3 


3. 


assign: 


V^n-/ C°^) = (°^) 



Xn-l = Xr> 

X n (‘^) - 



4 . 



E. 1. 



go to G.l. 



compute: 





0 =^) 



2, output: n(number of steps) 



Q (^) 

sea 

'^ext (oi') 
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.ttcrlng efficiency versus size parameter for sphere 
m = 1 . 54 ( 1-iO.OOl ) 



cu 

o 

CO 
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Size parameter 



cattering efficiency versus size paraxneter for slipere 

m = 1.54(1-10.001 ) 
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